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We study theoretically a quantum dot in the quantum Hall regime that is strongly coupled to a 
single lead via a point contact. We find that even when the transmission through the point contact 
is perfect, important features of the Coulomb blockade persist. In particular, the tunneling into the 
dot via a second weakly coupled lead is suppressed, and shows features which can be ascribed to 
elastic or inelastic cotunneling through the dot. When there is weak backscattering at the point 
contact, both the tunneling conductance and the differential capacitance are predicted to oscillate 
as a function of gate voltage. We point out that the dimensionless ratio £ between the fractional 
oscillations in G and C is an intrinsic property of the dot, which, in principle, can be measured. We 
compute £ within two models of electron-electron interactions. In addition, we discuss the role of 
additional channels. 



I. INTRODUCTION 

The Coulomb blockade occurs in an isolated meso- 
scopic island when the capacitive charging energy to add 
a single electron suppresses the discrete fluctuations in 
the island's charged In the past several years this physic* 
has been studied extensively in-both metallic systemscl 
and in semiconductor structures.u The Coulomb blockade 
is most easily probed by measuring transport through an 
island which is weakly coupled to two leads. By varying 
a gate voltage Vq, which controls the chemical potential 
of the island, peaks in the conductance are observed each 
time an additional electron is added. Between the peaks, 
the conductance is activated, reflecting the Coulomb bar- 
rier to changing the number of electrons on the dot. 

The coupling to the island via tunnel junctions intro- 
duces fluctuations on the dot which relax the discreteness 
of its charge. When the conductance a of the junctions 
is very small, a <C e 2 /h, this effect is weak and gives rise 
to the well known cotunneling effect, in which electrons 
may tunnel virtually through the Coulomb barrieriju 
Provided one is not too close to a degeneracy between 
different charge states, this physics may be described sat- 
isfactorily within low order perturbation theory in the 
tunnel coupling. 

For stronger coupling to the leads, the perturbative 
analysis is no longer adequate, and a quantitative de- 
scription of the problem becomes much more difficult. 
Based on general arguments, the Coulomb blockade is 
expected to be suppressed when a w e 2 /h,a since in that 
regime, the RC decay time of the island leads to an un- 
certainty in the energy which exceeds the Coulomb bar- 
rier. However, this argument is unable to predict quan- 
titatively the nature of the suppression of the Coulomb 
blockade. 

A particularly well suited system to study this physics 
is a semiconductor quantum dot at high magnetic fields. 
In the integer quantum Hall effect regime, the states near 
the Fermi energy of a quantum dot are edge states, and 



have a simple, well organized structure, which is insensi- 
tive to the complicating effects of impurities and chaotic 
electron trajectories. In this paper we shall study a quan- 
tum dot in the integer quantum Hall regime which is 
strongly coupled to a lead via a quantum point contact 
with one (or a few) nearly perfectly transmitting chan- 
nels, p 

In a recent paper, MatveevQ has considered a dot con- 
nected to a single lead by a nearly perfectly transmitting 
point contact. When the transmission of the point con- 
tact is perfect (zero backscattering), he showed that even 
in the presence of a substantial Coulomb energy U , the 
differential capacitance C = dQ/dVa is independent of 
Vq. Thus, the equilibrium charge Q on the island is max- 
imally un-quantized, since charge is added continuously 
as the gate voltage is changed. He then showed that 
the presence of weak backscattering at the point contact 
leads to weak oscillations in the differential capacitance 
as a function of gate voltage with a period corresponding 
to the addition of a single electron. These oscillations 
signal the onset of the quantization of the equilibrium 
charge on the island. 

In addition to the capacitance measurements,LfEI the 
fate of the Coulomb blockade in an island connected to 
a single lead via a point contact can be probed by trans- 
port, provided an additional lead is present. This leads us 
to the interesting possibility of simultaneously measuring 
the conductance (a transport property) and the capac- 
itance (an equilibrium property) of a quantum dot. In 
this paper, we build on Matveev's work and compute, in 
addition to the differential capacitance, the tunnel con- 
ductance and I-V characteristics for an island connected 
to one lead with a nearly perfectly transmitting point 
contact and connected very weakly to another lead. 

Specifically, we consider the quantum dot depicted 
schematically in Fig. [I]. The center region surrounded 
by the gates forms a quantum dot, which is connected 
to the leads on both sides. In the v = 1 quantum Hall 
regime, there is a single edge channel going around the 
dot. The contact to the left lead is a tunnel junction, 
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characterized by a small tunneling amplitude t, which is 
controlled by the voltage on gate A. The right contact 
is a nearly perfectly transmitting point contact charac- 
terized by a small backscattering amplitude v, which is 
controlled by the voltage on gate B. The backscattering 
at the point contact involves tunneling of electrons be- 
tween the opposite moving edge channels at x — L and 
x = —L, where x is a parameter specifying the spatial 
coordinate along the edge channel. 




FIG. 1. Schematic view of a quantum dot connected to 
leads on both sides. Negative voltage on the gates (shaded 
area) confines electrons in the dot. Thick solid lines show the 
edge channels formed by strong magnetic field B, where the 
direction of the electron motion is depicted by arrows. Dashed 
lines show tunneling paths, where the tunneling amplitude t 
and the backscattering amplitude v are controlled by voltage 
on gate A and gate B, respectively. 

A crude estimation of the Coulomb energy with the 
above geometry gives U = Ke 2 /eL, where K is a dimen- 
sionless geometrical factor and e is the dielectric con- 
stant of the semiconductor material. Another geometry- 
dependent energy scale is the level spacing between 
edge states in an isolated dot, AE — nhvp/L. While 
there are few reliable estimates of the edge state veloc- 
ity vf, the ratio AE/U is typically of the order 0.1 in 
GaAs/AlGaAs heterostructures in the integer quantum 
Hall effect regime. □ The smallness of this ratio is assumed 
throughout this paper. 

Though the equilibrium charge on the dot is not quan- 
tized when the transmission at the point contact is per- 
fect, we nonetheless find that important features of the 
Coulomb blockade remain in the tunneling characteris- 
tics. Remarkably, most of these features can be explained 
within the usual weak coupling model. Specifically, we 
find that at zero temperature, the Ohmic conductance is 
suppressed by a factor of (AE/U) 2 below its noninter- 
acting (i.e. U — 0) value. This suppression of the tun- 
nel conductance is precisely of the form predicted by the 



theory of elastic cotunneling through a one dimensional 
system. Evidently, the cotunneling theory is more gen- 
eral than its derivation within weak coupling perturba- 
tion theory suggests. In addition, when the temperature 
T or the voltage bias eV exceeds AE, we find behav- 
ior analogous to inelastic cotunneling. In particular, for 
AE <C T <C U the tunneling conductance is suppressed 
by (T/U) 2 . Moreover, for AE < eV < U, the tunneling 
current varies as V 3 . 

Inelastic cotunneling behavior has. also recently been 
discussed by Furusaki and MatveevL3 for a quantum dot 
strongly coupled to two leads. While similar in many re- 
spects, it is worthwhile to distinguish the present work 
from that in Ref. jl(|. In that reference, the two point 
contacts are considered to be independent. An electron 
passing through one contact can never coherently propa- 
gate to the second contact. Elastic cotunneling can there- 
fore not be described within that framework. In general, 
the elastic cotunneling will depend on the complicated 
transmission matrix of the sample. In contrast, in the 
quantum Hall regime, the edge channels have a particu- 
larly simple structure and directly connect the leads. If 
the phase coherence length can exceed the dot's dimen- 
sions, then elastic cotunneling should be present. This 
physics is correctly accounted for in our edge state model. 
In addition, the beautiful results in Ref. |10f| concerning 
the Coulomb blockade for spin degenerate systems is not 
applicable here, since the magnetic field destroys that 
degeneracy. 

In addition to the limiting behaviors discussed above, 
we find that at low temperature T <C AE, there can 
be nontrivial structure in the I-V characteristic for V ~ 
AE. In particular, we find steps in the differential con- 
ductance, dl/dV as a function of bias voltage. In the 
weak coupling limit ^-thc existence of such steps has been 
discussed by Glattlifl They are a consequence of inelastic 
cotunneling when the many body eigenstates of the quan- 
tum dot are discrete. Observation of such steps could 
provide a new spectroscopy of the low lying many body 
states of a quantum dot. Remarkably, as we shall show 
in Sec. Ill, these steps can remain sharp in the strong 
coupling limit, even though it becomes meaningless in 
that limit to speak of discrete single particle states. 

When the transmission through the point contact is 
less than perfect, charge quantization is introduced in the 
dot. The rigidity of the quantization grows with increas- 
ing amplitude of the backscattering at the contact. For 
weak backscattering, this is reflected in oscillations in the 
differential capacitance as a function of gate voltage, with 
an amplitude proportional to the backscattering matrix 
element v. We find that similar oscillations, proportional 
to v, should be present in the conductance. A compar- 
ison of these oscillations should therefore provide infor- 
mation about the intrinsic structure of the quantum dot. 
We focus on the ratio of the fractional oscillations in the 
conductance to that of the capacitance, 
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Go and Co are the average capacitance and conductance, 
whereas G\ and C\ are the amplitudes of the oscillations 
as a function of gate voltage. By considering the frac- 
tional oscillations, Gi/Gq and Ci/G), the dependence on 
the tunneling matrix element t necessary to measure G 
and the capacitive lever arm rj necessary to measure C is 
eliminated. Moreover, since in the weak backscattering 
limit both quantities are proportional to v, the depen- 
dence on v is eliminated by taking their ratio. Thus, £ 
measures an intrinsic property of the strongly coupled 
quantum dot which is independent of the details of the 
tunneling matrix elements. We have computed £ within 
a constant interaction model, in which the Coulomb in- 
teraction couples only to the total number of electrons 
on the dot. We find that f « 1.59 for T, AE < U, inde- 
pendent of T, AE and U. More generally, £ will depend 
on the specific form of the electron-electron interactions. 

With a few modifications, the above considerations can 
be extended to the case in which there are more than 
one well transmitted channels. We find that the low bias 
linear conductance is less suppressed as the number of 
channels AT increases. It has been recently shown that 
the suppression faotor becomes (AE/U) 2 ^ in a constant 
interaction modelLil However, if we take it into consid- 
erations that the interaction strength may differ within 
the same channel and between different channels, we find 
that the factor still has the form (AE/U) 2 . In addi- 
tion, the ratio £ depends sensitively on the form of the 
electron-electron interactions between different channels 
on the dot. It is equal to zero if different channels do 
not interact and grows with increasing strength of the 
inter-channel interactions. Thus, £ is a measure of the 
inter-channel interaction strength. 

This paper is organized as follows. In Sec. |TJ, we de- 
scribe the model and derive the I-V characteristic equa- 
tion. In Sec. Ill we consider a single channel system 



with a perfectly transmitting point contact. We compute 
the current and discuss various features of the result. In 



Sec. IV, we consider the effect of the weak backscattering 
to the conductance in connection with the effect to the 
capacitance. We compute £ within two specific models 
of electron-electron interactions. In Sec. [VJ, we generalize 
the results of Sec. Ill and Sec. IV to multiple channel 
systems. Finally, some concluding remarks are given in 
Sec. M. 



II. EDGE STATE MODEL 

We begin in this section by describing our edge state 
model of a quantum dot strongly coupled to a single lead. 
Here we will consider the case in which only a single chan- 
nel is coupled to the lead by a nearly perfectly transmit- 
ting point contact. Later, in section V we will consider 
the case where more than one channel is transmitted. 



In the absence of interactions, it is a simple matter 
to describe this system in terms of the free electron 
edge state eigenstates. Such a description is inconve- 
nient, however, for describing effects associated with the 
Coulomb blockade, which are due to the presence of a 
Coulomb interaction. The bosonization technique, how- 
ever, allows for an exact description of the low energy 
physics in the interacting problem.Ej 

At energies small compared to the bulk quantum Hall 
energy gap, the many body eigenstates are long wave- 
length edge magnetoplasmons. These may be described 
as fluctuations in the one dimensional edge density, n(x), 
where x is a coordinate along the edge. Following the 
usual bosonization procedure, we introduce a field <j){x), 
such that n(x) = d x (f>/2TT. The Hamiltonian, which de- 
scribes the compressibility of the edge may then be writ- 
ten, 



H = J dx^-(d x 



(2.1) 



The dynamics of the edge excitations follow from the Kac 
Moody commutation relations obeyed by <f>, 



Or 



2n 



= iS(x — x'). 



(2.2) 



Using ( |2.l| ) and ( |2.2| ) it may easily be seen that dtn = 
vpd x n, so that the edge excitations propagate in a single 
direction at velocity vp along the edge. In this language, 
the electron creation operator on the edge may be written 

as ^(x) = e l ^ x \ 



Equations (2.1) and (2.2) are an exact description of 
a single edge channel of noninteracting electrons. We 
may easily incorporate interactions into this description. 
Specifically, we consider a "constant interaction model" 
in which the Coulomb interaction couples to the total 
number of electrons N on the dot, which depends on the 
edge charge between x = — L and x = L, 
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[<P(L)-<P(-L)}. 



(2.3) 



The self capacitance and the coupling to a nearby gate 
can then be described by the Hamiltonian, 



H„ + Hr. 
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N 2 + enV G N, 



(2.4) 



where n is a "lever arm" associated with the capaci- 
tance coupling to the gate. Since the interaction is still 
quadratic in the boson fields, it may be treated exactly 
in this representation. 

Now we consider tunneling between the edge channels. 
We consider the left lead in Fig. [l] to be a Fermi liquid. 
Without loss of generality, we model it as another v = 1 
quantum Hall edge, ch arac terized by a boson field <f>i with 
a Hamiltonian Hq in (2.1). Tunneling from the left lead 
into the dot at x = is then described by the operator 



= e »[*(o)-cfc(o)] = 



(2.5) 
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where we define 9 = (f)(0) — (pi(0). For the reverse process, 
T^_ = exp — i9. The left point contact may be character- 
ized by a tunneling Hamiltonian, 



H f = t cos 9 



(2.6) 



where t is the tunneling matrix clement. 

Similarly tunneling at the right point contact involves 
transfer of electrons between x — —L and x = L. The 
Hamiltonian describing these processes can be expressed 



vcos[4>(L) — (j>(— L)] = vcos2ttN, 



(2.7) 



where v is the backscattering matrix element. 

It is c onv enient to eliminate the linear gate voltage 
term in ( |2.4| ) by the transformation N — > N — No, where 
iVo = er)Vc I (AE + U) is the optimal number of electrons 
on the dot. Our model Hamiltonian describing a quan- 
tum dot with a single channel coupled to a lead may then 
be written 



+t cos9 + v cos 2n(N - N ). 



(2.8) 



We will find it useful in our analysis to represent the 
partition function as an imaginary time path integral. 
The action corresponding to H is then given by 

Sq = — [ dxdr d x cj)(vFd x <fi + id T (f>). (2.9) 

47T J 

Since the remaining terms in the Hamiltonian depend 
only on 9 and N, it is useful to integrate out all of the 
other degrees of freedom. The resulting action, expressed 
in terms of 9(t) and N(t), is then given by 



+ J dr [t cos 9 + v cos 27r(iV - N )] (2.10) 

where w„ is a Matsubara frequency. G _1 is the inverse 
of the Green's function matrix and can be explicitly ex- 
pressed as 
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where AE = irvp/L. 

We now briefly develop the framework for our calcu- 
lation of the tunneling current. The I-V characteristic 
(or equivalently the tunneling density of states) may be 
computed using the action Stot in (2.10). Working per- 
turbatively in the tunneling matrix clement t, we may 
compute the tunneling current in the presence of a DC 
bias V using Fermi's golden rule. 
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where |m) is an eigen state of the unperturbed Hamil- 
tonian with energy E m . Since the sum on n is over a 
complete set of states, we can re-express the above equa- 
tion as 
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(2.13) 



(2.14) 



In order to compute P(E), it is useful to consider the 
imaginary time ordered Green's function, which may be 
readily computed using path integral techniques. 



T{t) = (T T e w ^e- l9 W) 



(2.15) 



The real time correlation function may then be deduced 
by analytic continuation. The two terms in the commu- 
tater in ( p. 14 ) lead to 



with 



P(E) = P > {E) - P < (E) 



P>'<(E) = J dte lEt T{r -^it± 



H 



(2.16) 



(2.17) 



Limiting behavior of the tunneling current may be de- 
duced analytically from the asymptotic behavior of V(t). 
In particular, at zero temperature, the Ohmic conduc- 
tance is proportional to the coefficient of the 1 /t 2 term. 
In addition, it is possible to compute P{E) numerically, 
as is described in the following section and in more detail 
in appendix A. 



III. POINT CONTACT WITH PERFECT 
TRANSMISSION 

In this section, we will consider the case where the 
transmission through the right point contact is perfect. 
In this limit, there is no quantization of the charge on 
the dot. Since charge may flow continuously through the 
point contact, there is no preferred integer value for the 
charge. This ma y be seen clearly from the Hamiltonian 
in equation (2.8), where for v = the dependence on TVo 
is absent. 

It follows that as the gate voltage is varied that there 
should be no oscillations in either the differential capac- 
itance or the conductance. However, we will show below 
that tunneling through a large barrier onto the dot is 



still blocked by the charging energy. This blockade is a 
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result of the fact that the tunneling electron has a dis- 
crete charge, which cannot immediately be screened. 

To compute the tunneling current, we evaluate Vq(t), 
where t he su bscript indicates that v = 0. Since the 
action (2.10) is quadratic, we may write, 



Po(r) 



exp 



(3.1) 
(3.2) 



wh ere G ee is the top left element of the matrix G defined 
in ( 2.1C| ) . This may be rewritten as 
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(3.3) 



where r c is the short time cutoff (which is of order the 
inverse of the cyclotron frequency) and 
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2ttU +U 1 -e 



(3.4) 



The first term in ( p.3[ ) describes the response for nonin- 
teracting electrons, U — 0. This gives us a purely Ohmic 
tunneling current / = Gjj=oV. Gu=a is related to the 
transmission probability of a free electron through the 
left barrier, Gu=o = {e 2 /h)T L , where T L oc t 2 . 

In the presence of interactions, the low bias linear con- 
ductance may be deduced from the long time behavior of 
Vq(t). In the strong interaction limit, U 3> AE, we may 
est ima te the limiting behavior of the exponential factor 
in ( |3.3| ) by noting that f(ui) is approximately a constant 
in each of the following three limits: 



if lo < AE, 
f{u) Pa { 1 if AE < w < U, 
ifZ7<w. 



(3.5) 



At zero temperature, we thus have, to logarithmic accu- 
racy, 



PoW « - ex P- 2 / — 
v t > Jae u 



(3.6) 



in the long time limit (r 3> AE^ 1 ). It then follows that 
the linear con ductance is suppressed. Performing the in- 
tegral in (|3.3|) exactly, we find 



g = Ci t Tl {-w 



(3.7) 



with ci pa 3.11. This should be compared with the theory 
of elastic cotunneling, which is derived in the case of weak 
tunneling through both barriers, Tl,Tr -C 1. When the 



dot is a one dimensional system (as it is for quantum Hall 
edge states), the result has been shown to be 



(3.8) 



Evidently, the AE/U suppression predicted in (3.7) re- 
mains valid all of the way up to Tr = 1. 

At finite temperatures or voltage s, A E <C eV, T <gi U, 
the lower limit of the integral in (3^6) is cut off by T 
and eV. The resulting tunneling current may thus be 
obtained by setting AE — and written, 



e 2 {eVf+A ^ T 2 

I = C2-rT L — V, 

h 



U 2 



(3.9) 



with c 2 = 2TT 2 e- 2 ^/3 pa 2.07 where C is Euler's con- 
stant. This result is, again, exactly in accordance with 
the theory of inelastic cotunneling, setting Tr = 1. 

An alternative interpretation of this suppression of the 
tunneling current has been pointed out in Ref. |l3f| . Sup- 
pose that an electron tunnels into the dot. The dot 
would minimize its electrostatic energy by discharging 
exactly one electron. According to Friedel sum rule, the 
number of added electrons, which is —1 in this case, is 
equal to 6/tt where 6 is the scattering phase shift of the 
one dimensionaLjchannel. As in Anderson orthogonal- 
ity catastrophe,E.a the suppression factor in the tunneling 
rate is related to the phase shift by 



dl_ 

dV 



oc e ' 



(3.10) 



where e = max(A£', T, eV) is the low energy cutoff and 



21'- 

7T 
-1] 
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Com bining ( [3.10] ) and ( [3.11 ) , we may reproduce ( 3/7 ) and 
( |3.9[ ). Note that ( |3.11 ) differs from the usual orthogonal- 
ity exponent, (6/ir) z , by a factor of two. This is because 
we are tunneling into the middle of a "chiral" system 
consisting only of right moving electrons (or equally the 
end of a one dimensional normal electron gas.) 

Finally, we note that in the high bias limit, eV 3> U, 
we recover the linear I-V characteristic with an offset 
characteristic of the Coulomb blockade, 



(3.12) 



This offset is a consequence of the fact that at short times, 
the electron which tunnels can not be effectively screened 
by the point contact. 

In addition to the limiting behaviors described above, 
we have computed the I-V characteristic numerically at 
zero temperature, as explained in appendix A. Fig. ^ 
shows the differential conductance dl/dV. What is most 
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striking is that there are sharp steps whose sizes are ap- 
proximately AE. For a nearly isolated dot, .-a similar 
phenomenon has been pointed out by Glattli,H and can 
be understood to be a consequence of inelastic cotunncl- 
ing through a dot with a discrete energy level spectrum. 



are the energies of the particle and the hole respectively. 
The electron eigenstates in the lead are labeled by k and 
k', and V\ k and Vi k > are matrix elements of the coupling 
Hamiltonian. If t\ — e 2 <C U/2, the above equation can 
be approximated 



32ir |VT ei - e 2 



1.0 



0.8 - 



0.4 



(b) - 


— AE=0 
.1U 
.3U 












// 
/'/ 

7 
/ 




/;/ " 
'/ 






j/i 

u 





0.0 



0.2 



0.4 



0.6 



1.0 



eV/U 



FIG. 2. Differential conductance dl/dV as a function 
of scaled bias voltage eV/U, computed in the absence of 
backscattering (Tr — 1) for various AE. Curves are scaled 
so that dl/dV = 1 in the high bias voltage limit. It is easy to 
see that the curves consist of steps of approximate size AE. 



An inelastic cotunneling process leaves a particle-hole 
excitation in the dot. As the bias voltage increases, the 
number of available particle-hole combinations also in- 
creases. Because of the discrete nature of the energy 
spectrum of the dot, this increase in number occurs dis- 
continuously at every AE/e of the bias voltage, which is 
manifested in the tunneling density of states or the dif- 
ferential conductance. However, the above explanations 
are not fully adequate in our model because the dot is 
strongly coupled to the lead. For perfect transmission, 
the line width of a single particle energy level is approxi- 
mately r ~ AE. It means that the levels are as broad as 
the level spacing and the steps are expected to be wiped 
out altogether. The reason for the apparent discrepancy 
is that the lifetime of a many body excited state (i.e. a 
particle-hole pair) can be much larger than the naive sin- 
gle particle lifetime. 

If the dot is weakly coupled to the lead, the lifetime of 
a particle- hole excitation may be easily calculated. First, 
we assume the excitation is relaxed only through the pro- 
cess in which both the particle and the hole tunnel o 
of the dot. In analogy with the theory of cotunneling,! 
we use Fermi's golden rule to estimate the decay rate, 
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where /i is the chemical potential of the lead and ei and e 2 
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where Ai ea( i is the level spacing of the lead and Tr ~ 
\2irV\ 2 1 AEA[ eai i is the transmission probability. We as- 
sume |Vife| 2 « |V2fe'| 2 ~ \V\ 2 is constant in the given 
range. 

In our model, a simple consideration of the I-V charac- 
teristic equation in (2.13) and flAl| ) shows that the width 
of the step risers is proportional to (AE/U) 2 . Since the 
width of the risers is directly proportional to th e line 
width of the energy levels, it is evident that Eq. (3.14) 
is valid even up to Tr = 1 with a possible numerical 
factor. We thus conclude that even in the presence of 
a perfectly transmitting contact, there can be long-lived 
excited states in the dot, whose decay is suppressed by 
the Coulomb blockade. 

So far we have assumed AE is constant. If we allow 
non-uniform level spacing, degeneracy in the particle-hole 
excitation energy is lifted and all steps but the first one 
split into several sub-steps. As the degree of degener- 
acy increases with the energy, more splittings occur at 
higher bias voltages, finally making it hard to distinguish 
between steps. 

Before we close this section, let us consider other re- 
laxation processes. It is only when all relaxation rates 
are less than AE that it is possible to experimentally 
observe the steps. This criterion is equivalent to saying 
that the inelastic scattering length, 1$, is much longer 
than the circumference of the dot. It is known that in- 
elastic scattering is strongly suppressed in the quantum 
Hall regimeEJ most likely due to the difficulty of conserv- 
ing both energy and momentum when scattering occurs 
in a one dimensional channel. These steps may thus be 
observable in the quantum Hall regime. 



IV. POINT CONTACT WITH WEAK 
BACKSCATTERING 

In this section, we consider the case where there is 
weak backscattering at the right point contact. As the 
contact is pinched off, fractional charge fluctuations in 
the dot are hampered and the discreteness of charge be- 
comes important. For weak backscattering v, there i s an 
energy cost, proportional to v in the Hamiltonian (2.7) 
for non- integral charge configurations. This gives rise to 
oscillations in physical quantities such as the capacitance 
and the conductance as a function of gate voltage. The 
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period of these oscillations corresponds to changing the 
optimal number of electrons on the dot Nq by one. 

For nearly perfect transmission through the point con- 
tact, we thus expect the conductance and the differential 
capacitance C — —edN/dVa to have the form, 



C= C + Cicos27riVo 
G = Go + d cos2ttN q , 



(4.1) 
(4.2) 



where the oscillatory components C\ and G\ are propor- 
tional to v. An intrinsic quantity, which is independent 
of the backscattering amplitude v, the tunneling matrix 
element t, and the capacitive lever arm r\ associated with 
the gate, is the ratio 



G\/Gq 
Cxi Co 



(4.3) 



Using the model we have developed so far, we can cal- 
culate the capacitance and the linear conductance per- 
turbatively in the backscattering matrix element v. The 
differential capacitance is given by C = TcP In Z / dV^ / -q , 
where Z is the partition function. As shown by Matveev, 
to leading order in v the differential capacitance has the 
form ( |4.1| ) with average and oscillatory components given 
by 



Co = 



AE + U 

G\ = vqe \ y ' la 



2ne 



AE + U 



(4.4) 
(4.5) 



The average (• • -) is with respect to the ground state of 
the unperturbed action Stot hi (2.10). 

In order to compute the conductance, we must calcu- 
late V{t) in the presence of the perturbation v cos 2ir(N— 
No). To the first order in v, we find that V(r ) = 
Vo (t) +V\ (r) cos 27riVo , where Vo (t) is given in (EO) and 



Vx{r) 



T (jlHT)-Ho)] q 



TTe mr)-0{0))) (S V ) Q 



(4.6) 



where S v — f a dr'v cos 2ttN(t'). This may be written as 

Vl {r)=vVo{T)e- 2 ^^ )2 ) a 

x J dr' {1 - cosh2vr [G 8N {t - r') - G 8N (t')}} , 

(4.7) 

where Gg n (t) is the off d iagonal element of the Green's 
function defin ed in (2.10), which may be computed ex- 
plicitly using ( 2.1l| ). 

In order to compute the lin ear conductance, we must 
compute the large r limit of (4/7). For Ut ^> 1 Gbn(t) 
decays as (Ut) -1 , so that the integral is independent of 
t. We thus find 



d = vGoe- 27T2 ( N{0)2 }a X 2 J dr' [1 - cosh27rG eA r(r')] , 

(4.8) 

where Go is the zeroth order linear conductance. 

Using (4.4), (4.5), and (p~g|), we obtain an exact ex- 



pression for the ratio £, 



AE + U 
At: 2 



2 J dr' [1 - cosh27rG eA r(r' 



)]. (4.9) 



In the limit T, AE U, the integral approaches a finite 
value which depends only on U. In this limit we find 



£ w 1.59. 



(4.10) 



The cancellation of input parameters like v, U, and AE 
may tempt us to suspect £ be a universal number being 
constant for all samples. As will be shown below, how- 
ever, Eq. (4.10) is true only within a constant interaction 
model, in which the dependence of the interaction on 
the spatial separation is ignored. In order to see how 
£ changes with different models, let us consider a more 
general model whose interaction Hamiltonian is given by 



1 

8t? 



dxdx' d x c/>(x)U(x, x')d x (j>{x'). (4.11) 



The constant interaction model is regained by assuming 
U(x,x') — U to be uniform. It is sufficient for our pur- 
pose to consider just another example. We can think of a 
local interaction U(x, x') = 2LUS(x — x')9(L—\x\), which 
is certainly an extreme limit to the other direction from 
the constant interaction model. The appropriate Green's 
function is given by 



Ge N (u! n ) = — ( 1 - e ) , 



and then using (4.9) we get 



(4.12) 



(4.13) 



Now it is clear that £ depends on the form of the electron- 
electron interaction. It is, however, noteworthy that the 
values of £ computed in two extreme limits are of the 
same order of magnitude. 



V. MULTIPLE CHANNEL SYSTEMS 

In this section we generalize the considerations in the 
previous sections to the systems of which right contact 
(nearly) perfectly transmits more than one channel. As 
the conductance of the contact increases with the num- 
ber of well transmitted channels Af, the shorter RC decay 
time allows a higher uncertainty in the energy of the is- 
land. Therefore it is natural to expect that the effect 
of the Coulomb blockade become weakened in multiple 
channel systems, which will be confirmed below. 
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It turns out that most of the qualitative considera- 
tions for the single channel systems can directly be ap- 
plied to multiple channel systems. Similar calculations 
as in Sec. [II show that in the absence of backscatter- 
ing, the low bias linear conductance is still suppressed 
below its noninteracting value, although the suppression 
is less strong if M is bigger. On the other hand, there are 
new features arising from the introduction of additional 
channels, on which we will focus in this section. 

For integer quantum Hall states with v > 1, the edge 
channels tend to be spatially separated. The tunneling 
will be dominated by the coupling to the nearest edge 
channel. As indicated in Fig. g this means that the tun- 
neling and backscattering will occur in different channels. 
Then the Hamiltonian for J\f channels is represented by 



i 

+ t cos 6>i + v cos 2?r (Njj - N ) , 



(5.1) 



where the notation is similar to that in Sec. Ill with the 
subscript denoting the channel number except that <pi is 
the boson field of the left lead. We have redefined No 
as the optimal number of electrons in channel Af alone, 
which in general depends on lever arms for all channels. 




FIG. 3. Schematic view of a quantum dot analogous to 
Fig. 1, with two well transmitted channels. The symbols are 
the same as in Fig. 1. Note that the tunneling through the 
left contact occurs in the outer channel whereas the backscat- 
tering at the right contact occurs in the inner channel. 

One of the simplest models to study a multiple channel 
system is a constant interaction model in which the in- 
teraction Hamiltonian depends only on the total charge 
in the island. The interaction may be explicitly written 



Hi, 



U 



5> 



(5.2) 



along the same lines as in S ec. III. In this case, the 
function f(u>) defined in (3^3) is given by 



U 1 



2tt\uu\+NU (l-e-^^j 

Note that the limiting values of f(oj) are 

f ifw<A£;, 
f(u) « I ^ if AE < u <C AfU, 
[0 ifAfU<u. 

It immediately follows that 



(5.3) 



dV 



6 '■ T 

oc —T L 



(5.4) 



(5.5) 



where e = max(A_E, eV, T). The exponent 2/Af has 
been de rivecL in some other pape rs in several different 
contexts.llirEj It is clear from (5.5) that the conductance 
is less suppressed if there are more channels. 

However, it turns out that the nonanalytic behavior 
with exponent 2/Af is correct only to the extent that the 
constant interaction model is valid. As explained below, 
this is due to a special symmetry of the charging energy 
with respect to redistribution of charge among the differ- 
ent channels. We consider an effective-capacitance model 
which is one step more general and has been introduced 
and developed by several authors to remedy some prob- 
lems with the constant interaction model.t£l This model 
assumes that the edge channels are capacitively coupled 
metal bodies and the Coulomb interaction energy de- 
pends on the number of electrons in each channel. The 
Coulomb interaction part of the new Hamiltonian can be 
written 



(5.6) 



where U is an Af x Af matrix which can be determined 
experimentally. In order to get a clear understanding of 
the effect of this generalization, let us consider a simple 
specific example of the electron-electron interaction, i.e., 



Ujj = u[AfaSij + (1 — a)]. 



(5.7) 



The diagonal component {Ma + 1 — a)u is the magni- 
tude of the interaction strength within each channel and 
the off-diagonal component (1 — a)u is that between dif- 
ferent channels. This matrix is chosen such that if the 
lever arms are all equal to unity, the total capacitance 



e ju in the limit 



The calculation of the tunneling conductance proceeds 



Ctot = dQ tot /dV G = e 2 . U.r J 
AE = 0, independent of a. Note that we regain a con- 
stant interaction model if a = 0. As a grows we move 
away from the model, finally reaching an independent 
channel model at a — 1, where different channels do not 
interact. As in (3.E), the limiting behavior of f(ui) de- 
fined in (|3.3[) is given by 
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if u < AE, 

i 1 if A.B < w < Afau, 

A«J » < Ma+l-a <<; w <<; ^ 

ifAAi<u>, 



(5.8) 



It is clear fro m the above equation that if a — ► (constant 
interaction), (5.4) is restored and we get the exponent 
2/ Af, as we discussed earlier. When a is not small, there 
is no appreciable range in which f(uj) = 1/A/", and the 
differential conductance has a different exponent, 2, i.e. 



dl e" fey 



(5.9) 



provided e = max(AS, eV, T) <C Afau. 

The above considerations can be generalized to a model 
with a generic matrix U. Even though f(u>) is a compli- 
cated function which depends on all matrix elements of 
U, there exists an energy scale u which corresponds to au 
in the above special model, such that 



f(ui) wl if AE < tJ < TV- 



(5.10) 



Then Eq. ( |5.9| ) is valid if £ <C A/w, except for the factor 
depending on a. In general, u = for a constant inter- 
action model and it measures how far the used model is 
awayixom the constant interaction model. Van der Vaart 
et alS3 have measured the matrix elements of U for two 
Landau levels confined in a quantum dot. Even though 
both contacts were nearly pinched off in the reference as 
opposed to our model, it is suggestive to estimate the 
magnitude of u using the experimental data. With their 
particular setup, they got Un = 800, U22 = 1175, and 
U12 = 650 (all in fieV). A simple estimation with these 
numbers gives Afu ~ 260[ieV 3> AE, T, eV, which sug- 
gests that ( |5.9| ) must be used rather than (5.5) in this 
case. It has to be admitted that this is a naive estima- 
tion considering the difference between their experimen- 
tal setup and our theoretical model. Opening up a point 
contact would in general reduce the strength of electron- 
electron interactioins in the dot and it would change 
the capacitance appreciably. However, even though the 
above estimation of u may be merely speculative at its 
best, we expect fl5.9| ) has to be true if different channels 
are weakly coupled, which seems more general than the 
constant interaction limit, in practice. 

It is now evident that the exponent is 2 /Af only for 
the constant interaction model. This is due to a special 
symmetry of the constant interact ion model, i.e., the in- 
teraction part of the Hamiltonian (5.1) is invariant under 
redistribution of total charge among the different chan- 
nels. The effect of the symmetry to the exponent can be 
most easily understood in te rms of Anderson orthogonal- 
ity catastrophe.E_l Eq. (3.10) can be directly used with an 
appropriately generalized definition of 7, i.e., 



(5.11) 



where S{ is the phase shift in channel i. It needs only 
a little consideration of electrostatics to figure out 
Following the argument in Sec. Ill, let us suppose that an 
electron has just tunneled into channel 1 through the left 
contact. The number of electrons discharged from each 
channel —Si/n depends on the form of the interaction, 
provided they satisfy the constraint YlSj/ir = — 1. If 
we work in a constant interaction model, because the 
Hamiltonian depends only on the total charge, from the 
symmetry Si/n = —l/Af for all i. Therefore 



Af' 



(5.12) 



On the other hand, if we use an effective-capacitance 
model and the system is safely away from the constant 
interaction limit (Afu ~S> AE), it is always energetically 
favorable to take a whole electron from channel 1 and 
have the exactly same ground state charge configuration 
as before. Then Si/n = — 1 and J,- = (i = 2 . . . Af) , so 
that 



7 = 2(-l) 2 = 2. 



(5.13) 



Then Eqs. (ph 
( |3~T0| ) and ( fUll ) 



and d5.5| ) are readily reproduced from 
The physical distinction between the 
energy scales u and u is thus clear. When an electron 
is added to the dot, (Afu)~ x , which corresponds to the 
RC decay time, sets the time scale for the total charge 
of the dot to return to its original value. However, even 
after the total charge has been screened, there may be 
some imbalance in the distribution of charge between the 
channels. (Afu) -1 sets the scale for the relaxation of this 
imbalance. In the constant interaction model, there is no 
Coulomb energy cost for such an imbalance, so u — > 0. 

Now let us consider the effect of weak backscattering. 
As in a single channel model, the introduction of weak 
backscattering v results in oscillations in the capacitance 
and the conductance. However, an important difference 
arises from the fact that the tunneling and backscattering 
occur in different channels. ■— ■ 
It has been— shown both theoreticallyll3 
and experimentally!!! that the period of the conductance 
and the capacitance oscillations increases with increasing 
number of well transmitted channels. This is because the 
oscillations arise only from the quantization of iV/v, the 
number of electrons in the backscattered channel. When 
there are many perfectly transmitting channels, many 
electrons must be added to the dot to increase Njy by 1. 

The analysis of the amplitude of the oscillations is a 
little more co mpl icated. We will again focus on the rati o 
£ defined in (4.3), using the model interaction in fl5.7| ). 
We assume AE, T -c Afau and all lever arms are taken 
to be unity. One may include the lever arms explicitly, 
but it does not change the result qualitatively. Along 
the same lines as in Sec. 



[V, the fractional capacitance 



oscillation may be written 
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Ci =vc -2^(NMof) (^\ 1 



Gn 



A/'' 



(5.14) 



The fractional conductance oscillation may also be writ- 
ten 

d = wG e~ 2 " 2 ( JVAf(0)2 )o x 2 JdT'[l- cosMttGw^t')] , 

(5.15) 

where the Green's function is given by 
G 0n (t) = (T T x (T)iV)v(O)) 
duj n e- lUnT 

(1 - a)usga{ui n ) 



(5.16) 



(27r|w n | + Mau)(2ir\u n \ + Mu)' 
We may easily compute Gon(t) in several limits, namely, 







G(w(t) 







■ Tr(l-a) 



if t = 0, 
if r — ► oo, 



(5.17) 



and it is monotonically interpolated in between. The 
above equation is not helpful if a ~ 1, but it is suf- 
ficient for our purpose which is to see how £ changes 
as the system moves away from the constant interaction 
limit. Since Gqn(t) measures the response of Nj^, a pe- 
riod of time r after an electron is added into channel 1, 
the physical interpretation of the above limiting behavior 
is clear. It takes a time period of order 2n/Mu for the 
total charge of the dot to return to its original value. 
Channel M contributes to this process by discharging 
(1 — a)/M of a n ele ctron, which can be read from the 
second line of (5.17). The reason it is proportional to 
1 — a is that the inter-channel interaction strength is pro- 
portional to 1 — a. After a time period of order 2ir/Mau, 
the charge in each channel returns to its original value, 
which is reflected in the vanishing G6w(t) in the lon g 
time limit. In order to estimate the integral in ( [5.15 ), 
we may make a crude approximation by substituting a 
square function for GgAr(r), i.e., GgAr(r) = — m(l — a)/M 
_a/ qm ? and Gejv( T ) — otherwise. Then we 

get 



Go 



ve 



2n 2 (NM0f} 8?r ( 1 I Q ) gin 2 7F ( 1 - g) 



Ma 



2M 



and finally 



27V 1 



■ sin 



Tr(l-a) 
2M ' 



(5.18) 



(5.19) 



This is a good approximation if a <C 1. This is a mono- 
tonically decreasing function of a and as is explained be- 
low, it is a consequence of the fact that the tunneling 
and the backscattering occur in different channels. A 
bigger a implies weaker inter-channel interactions and 



consequently a weaker effect of the backscattering to the 
conductance. At a = 1 (independent channels), we can- 
not use the above equation, but we know that the con- 
ductance oscillation would eventually vanish because the 
backscattering potential does not affect the conductance 
at all, and therefore £ = 0. At a — 0, (constant inter- 
action) Gi/Go diverges, and so does £. It is because 
Gi/Go diverges as (u/e)7T if the low energy cutoff e 
is small, suggesting that the perturbation theory break 
down. Without detailed calculations, one might have 
been able to infer it from the following physical argu- 
ment. In a constant interaction model, the total number 
of electrons in the dot J2i Mi is the only gapped mode 
and there are M — 1 combinations of M whose fluctua- 
tions are not bounded, leading to divergences in individ- 
ual terms in the perturbation expansion. Therefore, we 
need to sum up all higher order terms in order to obtain 
a correct result- In a series of recent papers, Matveev 
and FurusakiEHi3 have calculated both the conductance 
and the capacitance oscillations nonperturbatively in a 
spin-degenerate two-channel model, which they related 
to the multichannel Kondo problem. Their calculations 
show that the oscillations are no longer sinusoidal and 
the period becomes M times smaller so that the maxi- 
mum occurs each time an electron is added to the dot as 
a whole (not channel M alone). Such results, however, 
clearly apply only in the case where the degeneracy is 
guaranteed by a symmetry and hence should not apply 
in this quantum Hall system. 

Without qualitative changes, the above considerations 
can be generalized to an effective-capacitance model with 
a generic matrix U . As in the discussions of the differen- 
tial conductance dl/dV earlier in this section, an energy 
scale u, which plays the role of au, can be determined 
from the given matrix U. In most real situations of quan- 
tum Hall effect edge channels, £ is a finite quantity which 
can be numerically calculated in the effective-capacitance 
model if all matrix elements of U are known. 



VI. CONCLUSIONS 

In this paper we have shown that characteristics of the 
Coulomb blockade, which are normally associated with 
the weak coupling limit, persist to strong coupling to a 
lead via a single channel point contact. In particular, (i) 
we find the analogies of elastic and inelastic cotunneling 
in the {e/U) 2 suppression of the tunnel conductance, (ii) 
We find that particle-hole excitations on the dot can ac- 
quire a long lifetime due to a "Coulomb blockade" to re- 
laxation. This in principle could lead to observable steps 
in the low bias differential conductances as a function 
of bias voltage, (iii) The high bias behavior of the I-V 
characteristic has an offset, indicating the presence of a 
Coulomb gap. We find similar conclusions when multi- 
ple channels are transmitted through the contact, though 
the suppression of the Ohmic conductance is reduced. In 
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the special case of the constant interaction model, when 
there is no penalty towards redistribution of charge be- 
tween the channels, the exponent of the suppression is 
modified (e/J7)w. 

When the transmission through the point contact is 
less then perfect, the oscillations in the conductance and 
the capacitance may be characterized by the dimension- 
less ratio £. While £ is independent of the tunneling 
matrix elements, it depends on the precise form of the 
Coulomb interactions. For a single channel, we have 
computed it for two different forms of the interaction, 
and its value is of order unity. For multiple transmit- 
ted channels, its value depends even more sensitively on 
the inter-channel interactions, which is zero when differ- 
ent channels are independent, and grows with increasing 
strength of the inter-channel interactions. 
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APPENDIX A: INTEGRAL EQUATION 

FOR P > (w) 

In this-,appendix, we will use Minnhagen's integral 
equationuJ method to compute the spectral density func- 
tion Pq(E) for the perfect transmission case. Functions 
will be given the subscript to explicitly show that they 
are calculated in the absence of backscattering. The cal- 
culation of the imagi nary time Green function Po( r ) is 
straightforward from ( |3.2j ), and by analytically continu- 
ing it we get 



P o > (t)=Po(r^it + 0+) 



e v v_/ 'o exp / du 

/o w 



a(oj) 



(Al) 



where the average (■ • -) is evaluated over the unper- 
turbed action (t — v — 0). The function a(uj) is defined 



a(uj) 



UJ 

2> 

uj 



dte lut (0(t)0(O))„ 



<I^K)| 2 ) 

2 



(l^n)| 2 ) | 



7TUJ AE \ 



AE) 



~ 2^ sin £| (1- cos £§) 



(A2) 



Now we differentiate (Al) with respect to t and Fourier 



transform it. Then we finally get an integral equation 

ujP^iuj) = / dw' a{uj')P^(uj-uj'). (A3) 
Jo 
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